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We reconsider the effect of disorder on the properties of a superconductor characterized by a sign- 
changing order parameter appropriate for Fe-based materials. Within a simple two band model, 
we calculate simultaneously Tc, the change in residual resistivity Apo, and the zero-energy density 
of states, and show how these results change for various types of gap structures and assumptions 
regarding the impurity scattering. The rate of Tc suppression is shown to vary dramatically according 
to details of the impurity model considered. We search therefore for a practical, experimentally 
oriented signature of a gap of the s± type, and propose that observation of particular evolution of 
the penetration depth or NMR temperature dependence with disorder would suffice. 
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Introduction. Determining the symmetry and struc- 
ture of the superconducting order parameter in iron- 
based superconductors (FeSC) is one of the main chal- 
lenges in this new field [l| . The Fermi surface is usually 
given by two or three T = (0, 0)-centered hole pockets and 
two M= (tt, 7r)-centered electron pockets in the 2-Fe zone 
composed primarily of Fe 3d states. Repulsive intraor- 
bital interactions between hole and electron pockets lead- 
ing to spin fluctuations are often assumed to lead to a su- 
perconducting order parameter which changes sign over 
the Fermi surface to lower the overall Coulomb energy. 
The simplest version of this state, called the s± state, 
is described by an isotropic order parameter on each FS 
with opposite signs for electron and hole-like pockets ■ 
The state may be highly anisotropic and even exhibit 
gap nodes, but still be considered s± provided the av- 
erage sign on hole pockets is opposite that on electron 
pockets. On the other hand, other theories suggest that 
orbital fluctuations may dominate the pairing interac- 
tions in systems of this type, favoring a gap with equal 
sign on all pockets, denoted 5++ 0]. 

Surprisingly, it has proven rather difficult to defini- 
tively distinguish these types of gap structures exper- 
imentally, in part because phase-sensitive experiments 
are challenging due to surface properties; because of the 
multiband nature of the electronic structure; and because 
the s± and s++ "states" are symmetry equivalent, trans- 
forming both according to the Aig representation of the 
crystal point group. At this writing, three experiments 
offer indirect evidence in favor of the s± state: the nearly 
ubiquitous observation of neutron spin resonance features 
in inelastic neutron spectroscopy (INS) 043 , a quasipar- 
ticle interference scanning tunneling spectroscopy (STS) 
experiment in magnetic field [lol |. and a phase-sensitive 
experiment on a polycrystalline sample which relies on 
significant statistical analysis Kontani et al. have 



provided an alternate explanation for the neutron res- 
onance features within an scenario via a postulated 
energy dependence of the quasiparticle relaxation time. 
In addition, several references [12l - ll4| have called atten- 
tion to a "slow" decrease of Tc in chemical substitution 



experiments [14Ml7| . which is then ascribed to the nat- 
ural robustness against nonmagnetic disorder of an s+_|- 
superconductor . 

It is important to understand what is meant by "slow" 
and "fast" Tc suppression in this context. At one ex- 
treme we have situations in which Tc is not suppressed 
by nonmagnetic disorder at all. According to Anderson's 
theorem, the critical temperature of an isotropic conven- 
tional s-wave superconductor with a single band of elec- 
trons is unaffected by nonmagnetic scatterers. From this 
statement it follows immediately that the same occurs for 
two bands in an isotropic S-^.+ state (with equal gaps), 
but also in an s± state with no interband scattering. At 
the other extreme, we know that magnetic scatterers in 
a conventional isotropic superconductor su ppr ess Tc ac- 
cording to the Abrikosov-Gor'kov (AG) law [18[; it is well 
known that nonmagnetic scatterers suppress Tc at the 
same "fast" AG rate in a 2-band s± state, provided the 
two densities of states Na — Ni, and two gaps Aa = — At 
arc equal in magnitude, and the scattering is purely in- 
terband in nature. Any deviation from these assumptions 
will slow the Tc suppression rate relative to AG. There- 
fore between these two extremes lie many possibilities for 
Tc suppression behavior which depend on details of elec- 
tronic structure and the relative amplitudes of inter- and 
intra-band scattering. To say whether Tc suppression is 
"slow" or "fast" requires a more precise statement about 
the rate relative to observable measures of disorder, and 
compared to theoretical predictions. 

Several theoretical calculations of Tc suppression have 
discussed the pairbreaking effects of nonmagnetic scat- 
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terers on model multiband superconductors with general- 
ized s-wave order 12 , 28 1 . In fact the situation is gen- 
erally even more complicated than discussed above or in 
these works, since chemical impurities may do more than 
simply provide a scattering potential: they may dope the 
system, or alter the pairing interaction itself locally. We 
therefore believe (see also Refs. [l|) that measurements of 
Tc suppression relative to the amount of chemical disor- 
der are not particularly useful to determine gap structure 
in multiband systems unless a microscopic theory of the 
full influence of a single impurity can be developed, which 
seems challenging in the near future. Nevertheless, there 
are a few simple steps which can be taken to improve the 
situation. The first is experimental in nature: one needs 
to find a way to create pointlike potential scattering cen- 
ters, so as to create systematically disordered systems to 
which the above theoretical works apply. Based on expe- 
rience with many other materials, the closest approach 
to this ideal is achieved with low-energy electron irra- 
diation, which is thought to create interstitial-vacancy 
pairs. Experiments of this type are being performed cur- 
rently, and it is one of the goals of this work to make 
predictions to guide the analysis of such data. The other 
needed improvements are theoretical in nature: first, the 
pairbreaking theory must be extended to relate Tc only to 
directly measurable quantities, like the change in resid- 
ual [T — !■ 0) resistivity caused by the disorder, rather 
than to any theoretically meaningful but empirically in- 
accessible scattering rate parameter. Second, since the 
theory involves many parameters, the robustness of any 
claimed fit must be tested by the simultaneous prediction 
of other quantities which depend on disorder, such as the 
low-tcmperature penetration depth or NMR relaxation 
rate. Finally, it would be useful to have ab initio calcu- 
lations of vacancy and interstitial potentials to constrain 
the impurity parameters used. This has been attempted 

0, 



30 1 recently. 



for chemical substitucnts 

Model. We consider a system with two bands a and 
6. The ^-matrix equation in the two-band model has the 
form 



^ ^impt, 
t = U + Ugt, 



(la) 
(lb) 



where riimp is the concentration of impurities, t(ninip) = 

SLo*^*''®'^*' g('^imp) = go^'To + gi'X'n and ® represents 
a product of band (bold) and Nambu ( " ) matrices, go 
and gi arc local Green's functions in the tq and ri chan- 
nels (we have assumed particle-hole symmetry in order 
to neglect ga). Due to the translational invariance of the 
disorder-averaged system, g is diagonal in band space. 
We now assume a simple model for impurity scattering 
whereby electrons scatter within each band with ampli- 
tude V and between bands with amplitude u, 



u 



U V 



0T3. 



(2) 



The ^-matrix components are found from Eq. (jlb|) to be 



V 



gibu^ + giav'^ - gia (u^ - v'^Y Sg^ 



V 



(3) 



and =tS( a -n- b), where 



(4) 



V=l- {6gl + 5gl) + 5gl5gl 

-2v? {goagob - giagib) ■ 

with the abbreviation 5g^ = g 

Tc suppression. The linearized multiband gap equation 
near Tc is (see, e.g. Ref. |23 ) 



gL- 



A,(k) = 2T 



LOri—l^c 

E 

:',/3,tj„>0 



A,(kO 



Q/3 

kk' ~2 



(5) 



where ^/j is the linearized dispersion of band /3, and we 
introduced the shifted gaps and frequencies, A^(k') = 
A^(k') -I- Y^-p^ and = a;„ -|- iS^^-*. We will simplify the 
model above further in that we adopt a gap structure 
similar to that obtained from spin fluctuation theories: 
The gap on the (hole) pocket a is isotropic, A^, and 
the gap on the (electron) pocket b may be anisotropic, 
Afc = A" + /S.\{0), where 9 is the momentum angle 
around the b pocket and / dOAKO) = 0. The pairing 
potential is then taken as Vj^ = Vap(f>a{'i<.)4>i3{^'), with 



4>a — 1 + rSa.b cos 20, and 
tron pocket. 

This ansatz then gives three coupled gap equations for 

/^O Al^T 

equations in the compact form. 



is the angle around the elec- 



(Aa, A°, A^)"^ = A. In the A basis we can write the gap 



In 



1.13^ 



MA= CoMA, 



(6) 



the matrix A4_ = 
Co = ln(l.l3ff 



(1 + Vn V and the constant 

were introduced. Here V. is the inter- 
action matrix in the above basis. TZ is the orthogonal 
matrix which diagonalizes the matrix A, 

NbU^ -NbU^ 

-NaU^ NaU^ 

NbV^ + NaU^ 



A = 



V 



N 



(7) 



where Na is the density of states of the band a at the 
Fermi energy, 2?^ is Eq. (j4|) evaluated in the normal 
state, and ^ is a matrix with only diagonal elements. 



A, 



2ttTc 2ttTc 



A,; 



2ttTc 



where ^ is the digamma function and A^ being the eigen- 
values of matrix A. The maximum eigenvalue {Xmax{Tc)) 
of the matrix M_ is related to Tc, 

Tc ^I.UlUc c^p{- \ ). (9) 
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The sohition of Eq. ([9]) gives the value of Tc. 

Residual resistivity. The most direct observable mea- 
sure of scattering in suppression experiments is the 
residual resistivity change Apo, i-e. the change in the ex- 
trapolated T — !• value of the resistivity with disorder. 
We will assume that interference effects between elastic 
and inelastic processes are negligible, i.e. that the ef- 
fect on the p(T) curve when the system is disordered is 
essentially a T-independent shift upward. We therefore 
calculate Apo within the same framework as above, as- 
suming that all defects are pointlike. Our starting point 
is the result for the conductivity calculated from the so- 
lution of the Boltzmann equation or equivalently Kubo 
formula for isotropic pointlike impurity scattering with 
amplitudes given by Eq. ([2]) in band space. In the zero 
frequency limit (dc-conductivity), there are no interband 
transitions, and the total conductivity in a;-direction is 
the sum of the Drude conductivities of the two bands, 
a = (7a + <^b, with (Tq, = 2e'^Na{v'^ ,j.)Ta, where Va.x is 
the component of the Fermi velocity in cc-direction and 
Ta the corresponding single particle relaxation time ob- 
tained from the self-energy in t-matrix approximation 
Eq. dH), = -2ImI]L"^. Note that t'^ contains both 
contributions from the intraband and interband impu- 
rity scattering processes. The transport time and single- 
particle lifetime are identical within this model because 
of our assumption of pointlike s-wave scatterers, which 
implies that corrections to the current vertex vanish. A 
finite spatial range of the scatte ring potential will tend 



to steepen the Tc vs. Apo curve [31 



Results: vs. Apo- We now solve Eqs. ^ for T^, and 
calculate simultaneously the change in resistivity Ap^ at 
T — > 0. Unlike Tc vs. riimp or various scattering rates, 
Tc vs. Apo can then be compared directly to experi- 
ment. Clearly results will be parameter dependent, how- 
ever, so we here specify our precise assumptions regard- 
ing the electronic structure. For concreteness, we focus 
on the BaFe2As2 (Bal22) system on which the largest 
number of measurements have been reported. We as- 
sume a density of states on each Fermi surface sheet of 
Na = 3.6 and 2.7/unit ceh/eV/spin, for the "effec- 

tive" hole and electron pockets, respectively that approx- 
imately describes the imbalance in the densities of states 
that also has been seen with angle-resolved photoemis- 
sion (ARPES) [32-34| and is consistent with the density 
of states of Bal22 arising from Fe d-orbitals according to 
DFT calculations [sHi with an effective mass renormaliza- 
tion of z = 3. We take the root- mean-square Fermi veloc- 
ities as wf.o = 2/3x 10^ m/s and WF.b = 10^ m/s from Ref. 
[3^, Table I wp.i, and renormalize them by the same fac- 
tor of z = 3 to approximately match the velocities found 
in ARPES experiments |32[434l | . In the transport calcula- 
tion, the component of the Fermi velocities in the direc- 
tion of the current is taken to be ^) = 1/2 due 
to quasi-cylindrical Fermi surface. The pairing poten- 



tials are Ka = Vbb = 0.05, Vab = Vba = -0.04. For the 
isotropic case (?■ = 0) we obtain the zero temperature gap 
values of A°q = —1.79Tco and AJJq = 1.73Tco, whereas for 
the nodal case (r = 1.3) these are A°q = — 1.22Tco and 
A^Q = 1.23rco with the critical temperature chosen as 
Tco = 30K. We have fixed the intra-band scattering po- 
tential at an intermediate strength value oi v = 0.25, but 
show results for other values in the supplementary mate- 
rial. Potentials are given in eV and we set ft = fee = 1. 

In Fig. [1] we now exhibit Tc suppression vs. the cor- 
responding change in residual resistivity Apo as defined 
above, both for a fully isotropic s± gap (r = 0), as well as 
a gap which has nodes on the electron pockets (r = 1.3), 
for a range of ratios u/v. It is clear that a wide variety 
of initial slopes and critical resistivities Apg for which 
Tc — is possible, depending on the scattering charac- 
ter of the impurity. The variability of the suppression 
rate with the ratio of inter- to intraband scattering has 
been noted by various authors [i^,!!!] before this. In fact, 
Efremov et al. [1^ have shown that the various Tc sup- 
pression curves of the isotropic s± gap fall onto universal 
curves depending on whether the average pair coupling 
constant (A) <,=,> when plotted against the inter- 
band scattering rate. While these results are intriguing, 
they are difficult to compare directly with experimentally 
measurable quantities. Other works have made compar- 
isons with the resistivity changes (for example Ref. Il4l ). 
but have typically presented results for s± states only 
for a single set of impurity parameters corresponding to 
the fastest rate of Tc suppression. Such assumptions lead 
always to critical Apo's comparable to the smallest ones 
seen in Fig. [TJ of order lO's of p^lcm. Here we see that 
more general values of the parameters can easily lead to 
much slower Tc vs. Apo suppression rates by disorder, 
with critical disorder values of Ap§ of order mficm. As 
discussed in Li et al. [l4| , such values are typical of chem- 
ical substitutions on various different lattice sites; here we 
see that such "slow" Tc suppression does not rule out the 
s± state, even within the assumptions of our potential 
scattering only model. 

Results: density of states. Thus far we have focussed 
on Tc suppression itself, but it is clear that the parame- 
ter space is vast, so a real understanding of the effects of 
disorder in a given situation will probably depend on cor- 
relating the results of several experiments. Other quanti- 
ties which are quite sensitive to disorder are the tempera- 
ture dependence of the low-T London penetration depth 
AX(T) and nuclear magnetic spin-lattice relaxation time 
T^^. Within BCS theory, these quantities are controlled 
by the low-energy density of states. In the pure sys- 
tem, the nodal structure then determines the power law 
of temperature, and one generically expects AA(T) ~ T 
for gap line nodes except in very special situations [37l |. 
In the presence of a small amount of nonma gne tic dis- 
order, a finite density of states is created H, 3^ which 
leads automatically to a T^ term in the penetration depth 
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FIG. 1. (Color online) a) Normalized critical tempera- 
ture Tc/Tco vs. disorder-induced resistivity change Apo 
for isotropic s± wave paring for various values of the in- 
ter/intraband scattering ratio a = u/v. Insert: same quan- 
tity plotted over a larger Apo scale, b) Same as a) but for 
anisotropic (nodal) gap with anisotropy parameter r — 1.3. 



FIG. 2. (Color online) a) Schematic evolution of the density 
of states with increasing disorder for system with intra- and 
interband scattering, b) Top: Fermi level density of states 
Nb{0) (nodal band) as shown in Fig. [l]vs. Apo for various 
values of scattering ratio u/v. Bottom: for anisotropic s++ 
state with Vat identical in magnitude but > 0. 



371 140|. SuTiilar considerations lead to ^ in a 
pure nodal superconductors, but a linear T term induced 
by disorder. If the state is of s character, the gap nodes 
are not symmetry protected and can be lifted by further 
addition of disorder [24, 41|, thereby removing the terms 



in AA(r) or proportional to the disorder-induced 
density of states A^(0). 

In this work we note a further possibility in the dis- 
order evolution of the low energy DOS of a nodal multi- 
band s±-wave superconductor, namely that a reentrant 
behavior of A^(0) can occur after lifting of the nodes. 
The reason is that, in a situation dominated by intra- 
band scattering but with nonzero interband scattering, 
anisotropy of gaps on each individual sheet will be aver- 
aged by intraband disorder quickly. If the state is s± , a 
midgap impurity state can then be created by interband 
scattering, and grow until it overlaps the Fermi level, 
as shown schematically in Fig. [5] (a). The residual den- 
sity of states N{0) = -ImX;kTrG(k,w = 0)/(27r) (G is 
the Nambu Green's function) effectively determines low- 
energy thermodynamic behavior, so we have plotted it 
for the anisotropic band as a function of increasing dis- 
order in Fig. [21 for both s± and S-(-+ states. In the former 
case the reentrant behavior is clearly seen. 



The corresponding sequence in the s± penetration 
depth AA(r) would be T ^ ^ exp(-f7G/T) ^ T^, 
where Qq is the minimum gap in the system, while for 
the NMR spin- lattice relaxation rate Tj~^, the analogous 
evolution should he ^ T ^ exp(-rjG/r) T. In 
the s-(--|_ case, the last step in each sequence is entirely 
absent, since interband scattering cannot give rise to low- 
energy bound state formation. 

Realistic impurity potentials. It is clear from the above 
analysis that we have established that there is a wide 
range of possibilities for the behavior of T^, in an s± su- 
perconductor, as well as low-temperature properties like 
the penetration depth, when disorder is systematically in- 
creased. Within this range lie all current data on FcSC, 
from chemical substitution experiments which suppress 
Tc very slowly, to recent electron irradiation data which 
do so much faster. To make more precise statements, one 
needs to have some independent way to fix the scattering 
potential of a given impurity, and in particular the rel- 
ative proportion of inter- to intraband scattering. This 
has been attempted by a number of authors, with results 
which are not always in agreement, even qualitatively. 
In Kemper et al. [29[, the ratio between inter- and in- 
traband scattering was found to be of order a = 0.3 for 
Co in Bal22, which would lead according to Fig. [1] to a 
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critical resistivity strength of about 300 /ificm, roughly 
in accord with experiment 14 1. Onari and Kontani[13j 
have made the important point that the "natural" for- 
mulation for a model impurity potential, i.e. diagonal in 
the basis of the 5 Fe d-orbitals, automatically leads to 
significant interband scattering if one transforms back to 
the band basis. However, simple estimates show that de- 
pending on details a for onsite Fe substituents can vary 
between 0.2 and 1, again leading as seen in Fig. 1 to a 
wide variety of possible Tc suppression scenarios. 

Conclusions. We have argued that s± pairing cannot 
be ruled out simply because the suppression is "slow" 
according to some arbitrary criterion. The definitive ex- 
periments along these lines will most probably involve 
electron irradiation, where one can be reasonably sure 
that the defects created act only as potential scattcrcrs. 
In this case we find critical resistivities for the destruc- 
tion of superconductivity which vary over two orders of 
magnitude according to the ratio of interband to intra- 
band scattering. Results for the s± state are then not 
inconsistent with experimental data, but proof of sign 
change of the order parameter relies on knowledge of the 
impurity potential, which requires further ab initio cal- 
culations for each defect. As an alternative approach, 
we have proposed that systematic variation of disorder 
could give rise to a clear signature of s± pairing in the 
low energy Fermi level DOS N{0). In an s± state, A^(0) 
could increase with disorder, vanish again due to node 
lifting, and increase again after due to impurity bound 
state formation. This "reentrant" behavior of the DOS 
will be reflected in the temperature dependence of low- 
temperature quasiparticle properties like the penetration 
depth or nuclear spin relaxation time. For some ma- 
terials, this could be a "smoking" gun experiment for 
s± pairing. 
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[Supplementary Information] 



In this supplementary material wc give some details 
of how our results change when taking different values 
for the parameters impurity and pair potential parame- 
ters than uses in the main text, so that the reader may 
judge how robust our conclusions are. In the main arti- 
cle we have focused on the parent compound BaFe2As2 
(Bal22) and chosen values for the Fermi velocities and 
density of states at the Fermi level that are compatible 
with both DFT calculations [s^ and ARPES measure- 
ments (sil. [33! . However, there are still two parameters 
unfixed, namely the pairing potentials Vap and the im- 
purity potential v for scattering within bands (a full dis- 
cussion of the variation of the inter /intraband potential 
ratio a = u/v is included in the main text). Although 
the effective pairing potential Vap and average coupling 
constant_(A) = n^n^ Y.a,i3£{aM} NaVa/^Np as defined 
in Ref. |28| for our weak-coupling model, as well as the 
impurity scattering potentials u and v, are not known in 
experiments, our conclusions are consistent with different 
parameters within a reasonable range. If we increase v 
to V = 1.25 eV • Vc keeping all other parameters identical 
to those of Fig. 1 of the main text, the Tc suppression 
significantly slows, as seen in Fig. |3l with the exception 
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FIG. 3. (Color online) Tc/Tco vs. Apo for various val- 
ues of the inter/intraband scattering ratio a = u/v with 
V — 1.25 eV ■ Vc. a) for isotropic s± wave paring; b) for 
anisotropic (nodal) gap with anisotropy parameter r — 1.3. 



a=u/v 

FIG. 4. (Color online) The resistivity at half suppression 
a-s a function of a = u/v for various intraband impurity 
potentials v (measured in eV ■ Vc), other parameters taken as 
in the main article for the isotropic s± wave pairing (top) and 
for aniostropic (nodal) gap (bottom). 



of the value a = 1, which plays a special role in the 
theory of 2-band s± superconductivity, as can be eas- 
ily checked analytically. While Ref. y has argued that 
the interband scattering potential u should be generi- 
cally large for any chemical substituent, there is no rea- 
son to expect a = 1 to hold exactly, and therefore we see 
that large critical resisitvities ApQ arc even more likely 
to be found for stronger impurities (the unitarity limit 
w — ?> 00 with fixed a is pathological in this model[28| and 
we have not considered it here). The special role of the 
value a = I can be illustrated by plotting the resistivity 
Api/2 at which the critical temperature is suppressed by 
half, Tc = O.STco, as shown in Fig. 21 which may be com- 
pared with experiments. Note a ~ 1 yields the fastest Tc 
suppression independent of the impurity potential in the 
physical regime v > u. 

Finally, we also mention the effect of choosing other 
pairing potentials Va^ that lead to different values of 
(A). As explained in Ref. |28|, for isotropic s± paring, 
when Tc is plotted vs. the effective interband scattering 
rate, it follows three different universal curves accord- 
ing to whether (A) >,=,< 0. We have used a value 
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(A) = 0.037 in our paper. We have examined other 
parameter sets with negative (A), and found no essential 



difference in when plotted against the residual resis- 
tivity Apo, which of course depends on both intra- and 
interband scattering. 



